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We introduce a formalism of fractional diffusion on networks based on a fractional Laplacian 
matrix that can be constructed directly from the eigenvalues and eigenvectors of the Laplacian 
matrix. This fractional approach allows random walks with long-range dynamics providing a general 
framework for anomalous diffusion and navigation, and inducing dynamically the small-world 
property on any network. We obtained exact results for the stationary probability distribution, 
the average fractional return probability and a global time, showing that the efficiency to navigate 
the network is greater if we use a fractional random walk in comparison to a normal random walk. 
For the case of a ring, we obtain exact analytical results showing that the fractional transition and 
return probabilities follow a long-range power-law decay, leading to the emergence of Levy flights 
on networks. Our general fractional diffusion formalism applies to regular, random and complex 
networks and can be implemented from the spectral properties of the Laplacian matrix, providing 
an important tool to analyze anomalous diffusion on networks. 


PACS numbers: 89.75.Hc, 05.40.Fb, 02.50.-r, 05.60.Cd 


INTRODUCTION 


The recent burst of work around the idea of networks 
can be explained by the importance of this concept in 
a vast range of fields, which includes both the struc¬ 
tural features and the functional properties of networks 
[lJ-Q. In particular, we are interested in random walks 
taking place on networks, such as regular @-@|, random 
and complex networks 9|, and more recently temporal 
1C-I1, multiplex Ts)-ll[ , and interconnected networks 
16l I17j |. Random walks are useful to analyze problems of 
searching and navigability on networks, with applications 
to many different fields, such as the propagation of epi¬ 
demics, traffic flow, and rumor and information spreading 
0,1 ; for a recent survey, see Ref. p3]. 

In a recent paper, the usual navigation strategy of tran¬ 
sitions to nearest neighbors is generalized by allowing 
long-range navigation on complex networks using Levy 
random walks |l9[. This generalization allows transitions 
not only to nearest neighbors but to second-, third- or n- 
nearest neighbors. This new strategy was inspired by the 
study of Levy flights, where the lengths of random dis¬ 
placements obey asymptotically a power-law probability 
distribution 20|. These Levy flights generate anomalous 
diffusion 21] and have been used as searching and nav¬ 
igation strategies by animals {12-0] and also in human 
mobility and behavior 2^-37|. 

On the other hand, it is well known that one can study 
anomalous diffusion, and in particular Levy flights, using 
a fractional calculus approach BEHH. Likewise, we 
introduce in this paper a fractional approach applied di¬ 
rectly to the dynamics on networks to address the prob¬ 
lem of anomalous diffusion and long-range navigation. 
Our fractional analysis applied to general networks pro¬ 
vides us with a framework to deal with a richer dynamics 


on complex networks that includes, among other things, 
Levy flights fl9j |. 

This generalized long-range navigation can consider 
some common situations in real networks. For instance, 
in social networks one can use the knowledge of the net¬ 
work beyond our first direct contacts or acquaintances. 
Currently, using social networks sites, one can identify 



FIG. 1. (Color online) A random walker performing long- 
range transitions on a network can move directly from node 
i to node j with a transition probability given by ucfp\j (the 
index 7 refers to a fractional dynamics described in the text). 
The nodes are three degrees apart, that is, the geodesic dis¬ 
tance is three, as indicated by the dashed line. The geodesic 
distance is the integer number of steps of the shortest path 
connecting two nodes. Using this long-range dynamics one 
can contact directly long-distance nodes without the interven¬ 
tion of intermediate nodes and without altering the topology 
of the network. 
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the friends of your acquaintances (second-nearest neigh¬ 
bors) or the friends of the friends of your acquaintances 
(third-nearest neighbors) and so forth to search for in¬ 
formation, a job, an expert opinion, etc. In this fashion, 
one can contact people two or three degrees away of your 
friend directly, without the intervention of your friend. 
This situation, which we all use almost on a daily basis, 
corresponds to a long-range navigation on a network: a 
social network in this case. We illustrate this dynamics 
in Fig. [l] where we depict a general network of nodes 
and links among them. A random walker performing 
long-range transitions can move directly from node i to 
node j with a transition probability given by (the 

index 7 refers to a fractional dynamics described below). 
The nodes are three degrees apart, that is, the geodesic 
distance is three, as indicated by the dashed line. The 
geodesic distance is the integer number of steps of the 
shortest path connecting two nodes. In this long-range 
fractional dynamics on networks one can contact directly 
long-distance nodes without the intervention of interme¬ 
diate nodes and without altering the topology of the net¬ 
work. 

We start with an overview of the formalism associ¬ 
ated with diffusion processes and normal random walks 
on networks. Then, we extend the formalism to the case 
of fractional diffusion on networks and obtain a random 
walker defined by a transition probability matrix that al¬ 
lows a long-range dynamics. We deduce the stationary 
probability distribution of this process and, in order to 
study the efficiency of the random walker, we calculate 
the average return probability to a node and the aver¬ 
age global time associated to the capacity to explore the 
network. 


DYNAMICS ON NETWORKS 

We consider undirected connected networks with N 
nodes, described by the adjacency matrix A with ele¬ 
ments Aij = Aji = 1 (where i,j = 1,..., N) if there is 
a link between node i and node j, and A^ = Aji = 0 
otherwise. We take A it = 0 to avoid loops in the net¬ 
work. The degree of the node i is given by ki = ^u- 

The Laplacian matrix L is defined as Lij = Sijki — Aij, 
where Sij denotes the Kronecker delta; this matrix L is 
interpreted as a discrete version of the operator (—V 2 ) 
El]. Hence, the diffusion equation in a network takes the 
form [ 1 , 3 , 0 : 

L\m = -j t \m)- a) 

The vector \ip(t)) describes the system at time t, \ip(t)) = 
|ml, where {| 77 x)}JW 1 represents the canoni¬ 
cal base of [39|. On the other hand, random walks on 
networks are defined in terms of the modified Laplacian C 
with elements Cij = ^ = Sij — w^j , where w^j = 


are the elements of the transition matrix W of the nor¬ 
mal random walk on a network, describing transitions 
only to nearest neighbors with equal probability [&, [f|. 
For continuous time, the dynamics of the random walker 
is determined by the master equation Q: 

^ = ( 2 ) 

;=1 

where Pij(t ) denotes the probability to find the ran¬ 
dom walker in the node j at time t starting from the 
node i at t = 0. The master equation, Eq. ©, de¬ 
scribes a Markovian process with a stationary distribu¬ 
tion pj° = limt-^oo Pij(t), i.e. the probability for a walker 
to be in node j in the limit of large times. For a nor¬ 
mal random walk is given by p°° = „ N k] [3, Q ■ An- 

y m — 1 rCm 

other important quantity in the study of the diffusive 
transport is the average return probability defined by 
po(t) = jf 11 US 0 ]. From Eq. m it can be 

shown that p 0 (t) = Em=i expHCm*], where iCJilr 
are the eigenvalues of the modified Laplacian C [4| . 


FRACTIONAL DYNAMICS ON NETWORKS 


Having defined the Laplacian matrix L and the modi¬ 
fied Laplacian matrix C related to normal random walks, 
we introduce a generalization of these concepts in order 
to study the fractional diffusion on networks. For recent 
reviews of the fractional calculus approach to anomalous 
diffusion, see Refs, [gj I 20 I, l38|. 

We are interested in studying a generalization of Eq. © 
that reads: 


L 7 |M) = ~j t I i’it)) 


0 < 7 < 1, 


(3) 


where L 7 is the Laplacian matrix L to a power 7 , where 
7 is a real number (0 < 7 < 1 ). 

In the limit where 7 —> 1, we recover Eq. ©■ In the fol¬ 
lowing part we study the consequences of this definition 
and the characteristics of the random walks behind this 
dynamical process. 

Since L is a symmetric matrix, using the Gram- 
Schmidt orthonormalization of the eigenvectors of L, we 
obtain a set of eigenvectors {|\E'j)}^L 1 that satisfy the 
eigenvalue equation LI’Fj) = Pj\'&j) for j = 1 
and (’Iql’I'j) = Sij, where pj are the eigenvalues, which 
are real and nonnegative. For connected networks, the 
smallest eigenvalue pi = 0 and 0 < p m for m = 2,. .., TV 
41]. We define the orthonormal matrix Q with el¬ 
ements Qij = (il'Lj) and the diagonal matrix A = 
diag(0, p 2 ,..., pn)- These matrices satisfy L Q = Q A, 
therefore L = QAQ T , where Q T denotes the transpose 
of Q. Therefore, we have [42]: 


N 

L 7 = QA 7 Q T I'M ('IT 


(4) 
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where A 7 = diag(0, pj ,..., p? N ). In this way, Eq. © 
gives the spectral form of the fractional Laplacian matrix, 
and therefore, 


L 7 \^j) = I*] |*j) for j = 1,... ,N. (5) 


This result indicates that in order to treat the fractional 
dynamics we can simply calculate the spectrum {pj}jLi 
and then calculate {p~j}f =1 - 

On the other hand, in analogy with the matrix £, we 
define the modified fractional Laplacian matrix £^ with 
elements = (L 7 )y/(L 7 )j,. This modified fractional 
Laplacian is related to the dynamics of a random walker 
on a network determined by a fractional transition matrix 
W^ 7 ) with elements w^\j given by = 8ij — 

Therefore: 


w h) =6 - 
W j.(7) 


( 6 ) 


where we define the fractional degree of the node i as 
k = (L 7 )ij. Notice that w-E = 0. Also, the fractional 
transition matrix for 0 < 7 < 1 is a stochastic matrix 
that satisfies J2h=i w i^li = 1- On the other hand, from 
Eq. © in the case 7 = 1 , we recover the normal random 
walk with a transition matrix given by = Ay /fcj. 

Now, the corresponding fractional stationary distribu¬ 
tion p°° ( 7 ) of the random walker, from Eq. © , and using 
the detailed balance condition k < f l ^Pj°( 7) = k^pf° (7) 
0 , is given by 


Using Eq. © we obtain: 


N 


L 7 )i? = J^J 2 ^ 2 ~ 2COS ^ le 


7 „-i 4>i(i-j) 


( 8 ) 


1—2 

N 


= ^E( 2 _ 2 cos W 7cos ^( i_ ^ ( 9 ) 


i= 2 


with (f> m = 1). In Eq. ©, the imaginary part as¬ 

sociated with sin — j)] cancels out and only the real 
part remains in Eq. ©. Now we present this result in 
terms of the geodesic distance dij, defined as the integer 
number of steps of the shortest path connecting node i 
and node j. For the case of a ring: 


<f,, = r j| if 

’ I w - |i - j\ if 


\i-j\ = 0 , 1 ,..., LfJ , 


where [xj is the floor function that gives the largest in¬ 
teger not greater than x. The distances dij for a ring 
satisfy 


cos 




= COS 


2 7 r 


( 10 ) 


and using this result in Eq. we have for the fractional 
Laplacian of a ring: 


(L 7 )y = (2 - 2 cos I j>i ) 7 cos [<j>idij\ 


N 


N 


47 


1=2 

N 


= ^Vsin 
N 


27 


1= 2 


cos [4>id i: 


(ID 


On the other hand, using the fact that dij = 0 if i = j, 
we obtain from Eq. HID the fractional degree: 


pT{ 7 ) = 


(L 7 )i 


El= i(L 7 ), 


£ 1 ) 


Em=l k 


(7) ' 


( 7 ) 


This is a generalization of the result p°° oc kj for 
normal random walks discussed before, and is recov¬ 
ered from Eq. © when 7 = 1 . The general result that 
relates this stationary distribution with the mean first re¬ 
turn time still applies and reads: (tT 9 ) = l/p?°( 7 ) 0 . 


The modified Laplacian £ describes a random walker 
using a transition matrix that allows only the passage 
from a node to one of its neighbors (which we recover 
when 7 = 1 ). In what follows we show that in the frac¬ 
tional case, when 0 < 7 < 1 , the random walker can move 
using long-range navigation on the network. For this pur¬ 
pose, we start with the example of a ring (ID lattice with 
periodic boundary conditions), where the eigenvalues of 
the Laplacian matrix are p m = 2 — 2cos[27r(m — 1)/IV] 
and the corresponding eigenvectors have components 
(m|'I'/) = e - ' 7 ? 1 )( T " — - 1 )/-v/iV (where i = y/— 1) [Til . 


fc| 7) = (L 7 ) 


4 7 

N 


N 

y^sin 27 

1 = 2 



( 12 ) 


Now, introducing Eqs. CD and (| 12 |) in Eq. ©, we ob¬ 
tain the elements tc-E of the fractional transition matrix 
for a ring: 


w. 


(7) 

i-*j 


E^sin 27 

f(J-i) 

N 

COS 

[ 2 tt(Z- 1 ) , ] 
N ai J 

E^sin 27 

*( 1 - 1 ) 

N 



• ( 13 ) 


The result obtained in Eq. m reveals the nonlocal char¬ 
acter of the emergent process behind fractional diffusion 
on networks, where the transition probability depends 
explicitly on the global distance dij , at variance with the 
case of normal random walks, where the transition proba¬ 
bility allows transitions only to nearest neighbors. In the 
limit N 1, the sum in Eq. © can be approximated 
by an integral that takes the form: 


(L 7 k 


1 

2 tt 



( 2-2 cos eye idi ^dd, 


(14) 
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• Unvisited nodes 
O Visited nodes 

# Position at time t 


t=25 t=50 t=100 



1 


1 Random walk 
using regular 
dynamics 




Random walk 
using fractional 
dynamics 


FIG. 2. (Color online) Monte Carlo simulation of a discrete-time random walker on a network. We choose for clarity a tree, 
with N — 100 nodes, but the qualitative results stand for any network. The discrete time t denotes the number of steps of the 
random walker as it moves from one node to the next node on the network. This discrete random process is governed by a master 
equation with a transition probability matrix that gives the probability of moving from node to node. The dynamics starts at 
t = 0 from an arbitrary node. We show three discrete times t = 25, 50,100 for three values of the parameter 7 = 1 , 0.75, 0.5. 
Here, we depict one representative realization of a random walker as it navigates from one node to another randomly. The case 
7=1 corresponds to normal random walk leading to normal diffusion. In this case, the random walker can move only locally 
to nearest neighbors and, as can be seen in the figure, the walker revisits very frequently the same nodes and therefore the 
exploration of the network is redundant and not very efficient. The cases with 7 = 0.75, 0.5 correspond to a fractional random 
walk leading to anomalous diffusion. In this case, the random walker can navigate in a long-range fashion from one node to 
another arbitrarily distant node. This allows to explore more efficiently the network since the walker does not tend to revisit 
the same nodes; on the contrary, it tends to explore and navigates distant new regions each time. All this can be seen in the 
figure for different times, and allow us to make a comparison between a random walker using regular dynamics and fractional 
dynamics. 


which can be evaluated analytically and is given by (see 
Ref. [43[ for a discussion): 


that corresponds precisely to the Levy random walks 
introduced in Ref. (l9| . 


(L 7 ki = - 


ndij - 7)r(2 7 + 1) 

7iT(l + 7 + dij) 


sin(7T7), (15) 


where r(x) is the Gamma function. Using Eq. (11511 in 
Eq. © we obtain that for a ring: 


r 7 ) r(7 + i)r(dj 3 - -7) 

^ « r(-7)r(dy+7 + i)‘ 


(16) 


For i. j such that d \j 1, and using the asymptotic 

result r(n V a) ~ T(n)n“, for an integer n 1 , and a 
real ct, we arrive at the result: 


(7) 

n) 


,-(l+2 7 ) 

% 


(17) 


This asymptotic expression is not valid when 7—^0 
or 7 —> 1. To summarize, for a very large ring and 
very large geodesic distances between nodes, we have 
shown explicitly that the transition probability that 
emerges from the fractional dynamics is a power law 


This long-range navigation, based on power laws 
of the shortest paths introduced by Ref. [19| is valid 
for general networks and has been explored by other 
authors afterwards [13, 0 Si. For this long-range 
fractional dynamics, the transition probability matrix 
can connect arbitrarily distant nodes and, in this sense, 
the problem can be mathematically equivalent to an 
abstract complete weighted network. 

In order to illustrate the effect of the fractional dy¬ 
namics of a random walker on a network, in Fig. [2] we 
show the Monte Carlo simulation of a discrete-time ran¬ 
dom walker on a network. We choose for clarity a tree 
(network without loops Q) but the qualitative results 
stand for any network. The discrete time t denotes the 
number of steps of the random walker as it moves from 
one node to the next node on the network. This discrete 
random process is governed by a master equation with 
a transition probability matrix that gives the probabil- 
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ity of moving from node to node. Given the topology 
of the network, we calculate the adjacency matrix and 
the corresponding Laplacian matrix L of the network. 
Then we obtain its eigenvalues and eigenvectors that al¬ 
low us in turn to get the fractional Laplacian matrix L 7 . 
Finally, using Eq. ©> we determine the transition prob¬ 
abilities for different values of the parameter 7 . The dy¬ 
namics starts at t = 0 from an arbitrary node. We show 
three discrete times t = 25, 50,100 for three values of the 
parameter 7 = 1,0.75,0.5. Here, we depict one repre¬ 
sentative realization of a random walker as it navigates 
from one node to another randomly. The case 7=1 
corresponds to normal random walk leading to normal 
diffusion. In this case, the random walker can move only 
locally to nearest neighbors and, as can be seen in the 
figure, the walker revisits very frequently the same nodes 
and therefore the exploration of the network is redundant 
and not very efficient. The cases with 7 = 0.75,0.5 corre¬ 
spond to a fractional random walk leading to anomalous 
diffusion. In this case, the random walker can navigate in 
a long-range fashion from one node to another arbitrarily 
distant node. This allows us to explore more efficiently 
the network since the walker does not tend to revisit the 
same nodes; on the contrary, it tends to explore and nav¬ 
igates distant new regions each time. All this can be seen 
in the figure for different times, and allow us to make a 
comparison between a random walker using regular dy¬ 
namics and a fractional dynamics; see the Supplementary 
Material 461. 

Finally, it is important to stress that this general dy¬ 
namics given by the fractional Laplacian matrix L 7 has 
embedded the seeds of a long-range dynamics on net¬ 
works of any kind. Not only involving the shortest paths, 
but also trajectories of any kind in the network. That is, 
the fractional transition probabilities introduced in Eq. 
© contains a global dynamics in networks. Here lies its 
importance and the fruitful applications that it may have 
for many real processes in networks. 


bility, defined by p^\t) = ± E^Ii pWW, as: 


A 7 L 


1 N 


(19) 


where (j 7 -* = 0, 0 < (m ^ < 2 with m = 2,..., N are the 
real eigenvalues of /A 7 A Processes where p^it) decays 
rapidly to the lim^c*, (t) = TV -1 explore efficiently 
new sites of the network. 

To analyze in detail this average return probability, let 
us use as before the case of a ring. By using the spec¬ 
trum of the Laplacian matrix of a ring and the defini¬ 
tion of we can determine analytically the spectrum 
{C / <7 \}(=2 f° r the modified fractional Laplacian. With 
these eigenvalues and Eq. m, in the limit N — > 00, 
and approximating the sum by an integral as before, we 
obtain: 


*’<*>-s 


r 2i r 


3 —(2—2cos0 Ft/fcW 


< 10 . 


( 20 ) 


Here k ^ = — r( 7 ) r G 7 + 1 ) Sln ( 7r7 ') denotes the elements 

7rl (1+7) 

in the diagonal of L 7 given by Eq. (fl5l) . In the limit 
7=1 we recover the well-known analytical result for 
normal random walks where p^\t) = e _t /o(t), where 
I n {x) denotes the modified Bessel function of the first 
kind [47] . Not only for 7 = 1, but even for some ra¬ 
tional values of 7, we can obtain an analytical result 
for this quantity. For example, when 7 = 1/2 we have 
p ( o /2 \t) = Iq (^) — L 0 (y), where L„(a:) denotes the 
modified Struve function [48| . On the other hand, for 
t > 1, Eq. ® can be expressed as: 


r 

[-2 27 (9 27 1 

/ exp 
Jo 

k < 7 ) 1 _ 


d8, 


( 21 ) 


which can be evaluated analytically and takes the form 

0 : 


FRACTIONAL RETURN PROBABILITY 


pSr’m 


(/A 7 )) 27 
4-777 



( 22 ) 


Now, we analyze the continuous-time random walks 
defined by the fractional transition matrix WA). Using 
the fractional matrix /A 7 ) ; n the master equation ©, the 
fractional occupation probability p[^ (t) satisfies: 


dt 




1=1 


(18) 


We are interested in the efficiency of the dynamics, de¬ 
scribed by Eq. m, to explore the network. In the frac¬ 
tional case a treatment similar to the analysis of Eq. © 
allows us to obtain the average fractional return proba¬ 


where v = (2tt ) 2 ' y /and r(a, x) is the incomplete 
Gamma function [ 1481] . Therefore, the average return 
probability p^\t) for a ring, with N —> 00 in the limit 


t 1, decays as a power-law given by p[) 7l (t) ~ t 1 A 27 '). 
This result generalizes the case 7 = 1 for a normal ran¬ 
dom walk, where p^\t) ~ t” 1 / 2 [ 4 , 47]. 

Let us analyze now this quantity for general networks 
in more detail. In Fig. [3] we show the average frac¬ 
tional return probability p q 7 ' 1 (t) using the eigenvalues of 
the modified fractional Laplacian /A 7 ) anc j gq (| 1 9 [) . No¬ 
tice that the fractional dynamics explores more efficiently 
the networks in comparison with the normal case 7=1 
(i.e. Pq ' 1 {t) decays more rapidly to the asymptotic value 
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FIG. 3. (Color online) Average fractional return probability 
Po\t) as a function of time for different networks with N = 
5000 nodes and calculated using Eq. (HU. (a) A ring, (b) a 
tree, and (c) a scale-free (SF) network of the Barabasi-Albert 
type ^]. We used different values of 7 as indicated in the 
insets. 



0.0 0.2 0.4 0.6 0.8 1.0 
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FIG. 4. (Color online) Average time T vs 7 calculated from 
Eq. (1231) for different types of networks with N = 5000 nodes, 
(a) Large-world networks: ring (ID lattice), the regular net¬ 
work with degree k = 4 obtained by the Watts-Strogatz (WS) 
model with rewiring probability p = 0 [stil l, a tree, a 2D lat¬ 
tice with dimension 50 x 100 and periodic boundary condi¬ 
tions. (b) Small-world networks: the WS network used in (a) 
with rewiring probability p = 0.7, a scale-free (SF) network of 
the Barabasi-Albert type [49}, and Erdos-Renyi (ER) network 
with p = log N/N [5l| . 


N 1 ). In particular, in Fig. [Ua) we observe the power- 
law decay of p^ oc f -1 /! 2 '*'), as predicted by our analyti¬ 
cal results for a ring, before reaching the value iV -1 . The 
effect of the efficiency of navigation due to fractional dy¬ 
namics is more noticeable for large-world networks (ring 
and tree) than for small-world (scale-free) networks. 

In order to quantify the efficiency to explore the net¬ 
work, we introduce a global time T as: 

,00 1 N 1 

T= (Po 7) (t)-Po 7) (oo))dt = — (23) 

Jo 3 =2 C j 

The inverse of the eigenvalues (j 7 '* are the characteristic 

times Tj 7 ' 1 that dominates the dynamics; these times are 
also relevant for the problem of synchronization in net¬ 
works [ 3 }. Thus, this global time T is the average of these 
characteristic times in pg 7 ' 1 (t ). In the context of Marko¬ 
vian processes r = NT is the Kemeny’s constant, that 
for random walks is the global time r = ^2jLiPj°{Tij), 
where (Tj) is the mean first passage time (MFPT) de¬ 
fined as the mean number of steps taken to reach the node 
j for the first time starting from the node i @, [HI |H2, j53j. 

In Fig. [4] we show the global time T as a function 
of 7 for different types of networks. In Fig. Eta), for 
large-world networks, the effect of the fractional dynam¬ 
ics reduces several orders of magnitude the time T, in 
comparison with the case 7 = 1 . In Fig. Et^) we show 
that even for small-world networks the fractional dynam¬ 
ics improves the efficiency to explore the networks. 


CONCLUSIONS 

In summary, we have introduced a formalism of frac¬ 
tional diffusion on networks based on a fractional Lapla- 
cian matrix that can be constructed from the spectra and 
eigenvectors of the Laplacian matrix. This fractional ap¬ 
proach allows random walks with long-range dynamics 
providing a general framework for anomalous diffusion 
and navigation in networks. We obtained exact results 
for the stationary probability distribution, the average 
fractional return probability and a global time. Based on 
these quantities, we show that the efficiency to navigate 
the network is greater if we use a fractional random walk, 
in comparison to a normal random walk. For the case of 
a ring, we obtain exact analytical results showing that 
the fractional transition and return probabilities follow a 
long-range power-law decay and, thus, the emergence of 
Levy flights on networks. Our general fractional diffusion 
formalism applies to regular, random, and complex net¬ 
works and can be implemented from the spectral prop¬ 
erties of the Laplacian matrix, providing an important 
tool to analyze anomalous diffusion on networks. Our 
results show how the long-range displacements improve 













































7 


the efficiency to reach any node of the network inducing 
dynamically the small-world property in any structure. 
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